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Abstract. In this talk, we shortly report results from our recent global DGLAP analysis of 
nuclear parton distributions. This is an extension of our former EKS98-analysis improved with 
an automated \ 2 minimization procedure and uncertainty estimates. Although our new analysis 
show no significant deviation from EKS98, a sign of a significantly stronger gluon shadowing 
could be seen in the RHIC BRAHMS data. 



1. Introduction 

The global analysis of nuclear parton distributions (nPDFs) is driven by the experimental fact 
that the deep inelastic structure functions F2(x,Q 2 ) measured from nuclear targets show a 
significant deviation from the free proton ones [TJ [2] . 

Perhaps the most simple theoretical approach to this observation is to make use of the 
factorization theorem of QCD that has proven to provide excellent description of inclusive cross- 
sections in free-nucleon collisions. In this approach, the cross-sections are of the generic form 

*AB->h+x = £ ffouQ 2 ) ® fj B &, Q 2 ) ® a l+ ^ h+x , (1) 

where where a - i +j—*h+ x [ s the perturbative QCD (pQCD) matrix element squared and /jS are 
the non-perturbative parton densities whose scale evolution obeys the DGLAP equations [3J. 

The purpose of the global analysis of nPDFs is to find out whether the observed differences 
in the structure functions, the nuclear modifications, can consistently be absorbed in to the 
input parton densities — in other words, do the nuclear modifications effectively factorize. If 
they do, the resulting nPDFs are of great practical interest, since they can be used as a input 
in any process that can be factorized as in eq. (JT]). In this framework the deep question about 
the dynamical origin of nuclear modifications is not addressed, on the contrary, one must be as 
unbiased to any model as possible. 

Indeed, three independent groups have shown that this approach works quite well: 

• EKS98 [U [5] was the first global analysis demonstrating that using the leading-order (LO) 
pQCD formalism, requiring momentum and baryon number conservation, one can reproduce 
the data from measurements of deep inelastic lepton-nucleus scattering (DIS) and Drell-Yan 
dilepton production (DY) in proton-nucleus collisions. The fit was done only by eye. 
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• HKM [6] & HKN [7J which are LO QCD analyses as well, were the first ones to exploit the 
automated ^-minimization procedure and to estimate nPDF uncertainties. 

• nDS [8] brought the global analysis of nPDFs to the next-to-leading order level in pQCD. 

Based on our earlier work, EKS98, we have performed a global reanalysis of nPDFs published 
in [9] and reported here. What was new compared to the EKS98 were the automated x 2 - 
minimization procedure and a first attemp for estimating the uncertainties. Although our main 
objective was to see whether we can improve the EKS98-fit and study the uncertainties, this 
was anyway a necessary stepping stone for us before we can extend our analysis to NLO-level. 

Motivated by the BRAHMS data [10| on inclusive hadron production in D+Au collision, 
which show a systematic suppression relative to p+p at forward rapidities, we raise an intriguing 
question about the possibility of having clearly stronger gluon shadowing than what has been 
hitherto seen in the global DGLAP analyses. 



2. The framework 

We define the PDFs fi(x, Q 2 ) of bound protons in a nucleus with mass number A as 

f l A (x,Q 2 ) = Rf(x,Q 2 )f^ L1 (x,Q 2 ), (2) 

where j^ CTE Q 6L1 re fers to the latest free proton PDFs by the CTEQ collaboration For 
the bound neutrons we assume the isospin symmetry d pr oton = ""neutron and vice versa. 
What we actually parametrize, are the nuclear modifications Rf(x,Q 2 ) at the initial scale 
Ql = 1.69 GeV 2 . At present, the lack of precision data forces us to consider only three different 
modifications: Ry(x, Qq) for all valence quarks, R A (x, Qq) for all sea quarks, and Rq(x, Qq) for 
gluons. 

The fit functions are parametrized in three pieces (c.f. Fig. [3]): 

R A {x) = c A + (c? + 4x)[exp(-x/x£) - exp(-x^/xf)], x<x A 

R A {x) = a A + afx + a A x 2 + a A x 3 , x A < x < x A (3) 

r>A/ \ ^0 — b A X j{ ~ 

^3 \ x ) = J: 7^4 > x e — x i 

(1 — X)P 

The first one covers the region from shadowing to anti-shadowing maximum at x„ , the second 
comes down to EMC-minimum at x A , and the third is for the Fermi- motion part. 
The ^4-dependence of the fit parameters is assumed to follow a power law 

z ? = z t Tet {^Y H , (4) 

where we have chosen Carbon (^4 re f = 12) as a reference nucleus. 

After fixing the continuity of Rfs and their first derivatives at if and x£, one is still left with 
42 free parameters of which the baryon number and momentum conservation eat only 4 away. 
This was still too much in order to obtain converging fits, and lots of manual work was needed 
too see which parameters were the most relevant ones. At the end there was 16 fit parameters. 



3. Results & Error analysis 

The experimental input in our analysis was about 500 points of DIS- and DY-data in a form 
Tida tu /dQ 2 dx \d<j P p Y / dx 2 dQ 2 



These covered 11 different nuclei from Helium up to Lead. Figs. [T] and [2] show some of these 
data and comparison to our result from minimization of 
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Figure 1. Calculated R$ 2 (x,Q 2 ) (filled symbols) are compared to SLAC [14J, E665 [15], NMC 
95 p2] and reanalysed NMC 95 data [H]. The asterisks denote our results calculated at the 
initial scale Qq, these are for the smallest-x data points whose scales lie in the region Q 2 < Qq- 



The number that characterizes the goodness of the fit is x 2 /d-o.f. (d.o.f = iVdata — 
parameters) which should be less than 1 if the fit is any good. In our case x /d.o.f. ~ 0.8 
indicating that the theory fits the data very well and that there is no serious sign why 
pQCD could not be trusted — within the considered kinematical range — also in the nuclear 
environment. 

The obtained nuclear modifications at the initial scale for Lead are shown in Fig. [3l This 
figure also presents our uncertainty estimates based on Hessian method of quantifying the 
uncertainties [12J, in which one expands the x 2 around the minimum w.r.t fit parameters £ 
as 

Ax 2 = X\i + *0 " X\0 = E (6) 
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Figure 2. Calculated i?^ y (x,Q 2 ) and F 2 Sn /F 2 c (filled symbols) are compared to E772 [13] and 
NMC data UE\ 



and the uncertainty of any quantity F(£) depending on PDFs is then obtained from 

[<5F(|)] 2 = A X 2 £ 



-i(9F(j)\ 
\ ^ ) 



(7) 



For an ideal x 2_ distribution Ax 2 ~ 18, which we used too, corresponds to "one sigma"-error. 

Due to technical difficulties obtaining converging fits, the EMC minimum of gluons and sea 
quarks was fixed to follow the valence quarks at Qq. This resulted as an unreliably small error 
bands for gluons and sea quarks, and they had to be computed separately. They are the "Large- 
x errors" in Fig. [3l The combined uncertainties are shown as a yellow bands. Interestingly, the 
old EKS98 parametrization lies within these uncertainties and there is no reason to release a 
new parametrization — EKS98 works just fine. 

However, one should be very cautious about these error bands! First, below x ~ 10~ 2 
there is no experimental data above Q 2 = 1.69 GeV 2 , and the behaviour at small-x region is 
constrained only by the sum rules and is bound to the form of the fit function. Second, the PDFs 
themselves depend on choices and conventions, like kinematical cuts, choosing the factorization 
scale, treatment of heavy quarks, and choosing and weighting data sets in forming x 2 - How these 
choices affect the obtained PDFs is not seen in the error analysis performed here — the error 
bands only reflect the experimental data and their errors. Third, there is no well-established 
way to choose Ax 2 , and our choice Ax 2 ~ 18 is probably quite restrictive. 




Figure 3. Initial nuclear ratios for Lead together with their uncertainty estimates 



4. Are the gluons different? 

The DY and DIS data leave the gluons very unconstrained — gluon dependence comes only 
through the DGLAP-evolution. One possible way to constrain the gluon sector could be the 
inclusive hadron production at RHIC. Figure [U shows data from the BRAHMS collaboration 
[TO] for -RdAui the ratio between charged hadron production in D+Au and p+p collisions as a 
function of hadron's pr- From the point-of-view of pQCD, this is computed via 
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Figure 4. Minimum bias inclusive hadron production cross sections in D+Au collisions divided 
by that in p+p collisions at \/sjV7V = 200 GeV at RHIC. The ratio -RdAu is shown as a function 
of hadrons transverse momentum at four different pseudorapidities. The BRAHMS data [TO] 
are shown with the statistical error bars and the shaded systematic error limits. A pQCD 
calculation for h + + h~ production with the nuclear modifications from present work and KKP 
fragmentation functions is shown by the black lines, and that with the strong gluon shadowing, 
shown at right, by green lines. 



E fi(?uQ) ® // (*2, Q) ® a l+ ^ k+l ® D k ^ h+X (z, Qf), (8) 

ijkl 



where the D(z,Q 2 )s are the fragmentation functions. In Fig. 0] the pQCD calculation which 
employs the nPDFs from the present analysis and KKP [19] fragmentation functions, is shown 
by dark line. It is evident that, especially at very forward direction, our prediction is above the 
data and the shape is not well reproduced. 

Since hadron production at the kinematical corner of Iow-pt and forward rapidity reaches 
the small- x region of PDFs where the gluon distributions are the dominant ones, this could be 
signaling a larger uncertainty in our gluon modifications than seen in Fig. O For this reason 
we present an example of gluon shadowing that is much stronger. This is shown on the right 
hand side of Fig. 01 and it really helps: using this gluon modification the corresponding curve 
for -RdAu is brought clearly closer to the BRAHMS data. 

However, it's still too early to draw very strong conclusions from this observation, but a 
systematic study in the context of global DGLAP analysis is needed |2L)j . One should also bear 
in mind that here we are considering effects at very low-pT region, where the simple LO pQCD 
picture is pushed to its very limits — although one can argue that in ratios like RdAu some higher 
order effects would partially cancel. Anyway, one should be very careful when interpreting these 
results. For example, it has been conjectured, that this particular BRAHMS data set could be 
a sign of parton saturation at work |21j . 
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